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0 . ABSTRACT 


A cyclic  queueing  network  with  two  servers  and  a finite 
number  of  customers  is  studied.  The  service  times  for  server  1 
form  an  EARMA(1,1)  process  (exponential  mixed  autoregressive 
moving  average  process  both  of  order  1)  which  is  a sequence 
of  positively  correlated  exponential  random  variables;  the 
process  in  general  is  not  Markovian.  The  service  times  for 
the  other  server  are  Independent  with  a common  exponential 
distribution.  Limiting  results  for  the  number  of  customers 
in  queue  emd  the  virtual  waiting  time  at  server  1 are  obtained. 
Comparisons  are  made  with  the  case  of  independent  exponential 
service  times  for  server  1. 
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1.  INTRODUCTION 


Relatively  little  is  known  about  queues  for  %^ioh 
interarrival  tines  and/or  service  times  are  not  independent 7 
some  studies  are  available  however  (cf.  ^inlar  [1967], 

Loynes  [1962],  Pearce  [1967],  Purdue [1975] ) . 

One  reason  for  this  seems  to  be  the  lack  of  tractable  models 

for  dependent  sequences  of  random  variables.  Recently,  models 
have  been  developed  for  sequences  of  dependent  exponential 
random  variables  (cf.  Jacobs  and  Lewis  [1977]).  These  models 
are  parametrically  relatively  simple  and  are  in  general  not 
Markovian.  Fortunately,  it  seoui  to  be  easy  to  model  various 
types  of  dependence  in  queues  using  these  sequences. 

The  dependent  sequence  of  exponential  random  variables 

is  defined  as  follows.  Let  {e^>  be  a sequence  of  independent 

randkxi  variables  each  with  an  exponential  distribution  with 

■wan  1*^,  0 < 1 < *.  Let  {J^}  and  {K^}  be  independent 

sequences  of  independent  {O,l}>random  variables  such  that 

P{J  > 1)  > 1-6  and  P{K  >1)  ■ l>p  where  0 < $ < 1 and 
n 11  * — 

0 ^ p < 1 are  fixed  constants.  For  n ■ 1,2,...  put 

<»■»>  *»-•*.♦  'nVl  ' 

idiere 


(1.2) 


♦ V. 


Tlw  enqueues  (X.)  le  eelled  an  BhMhd,!)  process  (ei^onential 

n 

mimed  mevinq  ewereqe  evtmreqeeeive  both  of  order  1)  and 
ie  eeltdd  tihe  antoreereelve  pert  of  x.. 


We  will  always  assume  that  has  an  exponential 

distribution  with  mean  independent  of  all  the  other  random 

variables.  Under  this  assumption,  Jacobs  and  Lewis  [1977]  have 
shown  that  (X^)  is  a stationary  sequence  with  exponential 
marginal  distribution  having  mean  A**^.  Furthermore,  the  X^'s 
are  positively  correlated  with 

corr(X^,X^^^)  - p^’^(l-e) [e(l-p)  + (l-e)pl  , k > 1 . 


If  6 > 1,  or  6 > 0 and  p « 0,  then  {X  } is  a 

n 

sequence  of  independent  randosi  variables.  The  process  {X^} 
is  in  general  not  Markovian  although  it  is  if  3 » 0,  in  which 
case  it  is  called  an  BAR(l)  process  (exponential  autoregressive 
of  order  1) . 

In  this  paper  we  will  consider  the  simple  case  of  a 
closed  queueing  network  with  two  servers  and  a fixed  number  of 
cycling  customers  N.  It  is  described  as  follows:  when  a customer 
finiidies  service  at  server  1 he  joins  the  end  of  the  queue  at 
server  2;  when  he  finishes  service  at  server  2,  he  rejoins  the 
end  of  the  queue  at  server  1;  the  service  discipline  is  first- 
in-first-out . 

This  closed  queueing  network  has  been  used  in  computer 
studies  to  model  multiprogranmed  ocMuputer  systems  (cf.  Gaver 
and  Shadier  (1971))  and  one  is  interested  in  obtaining*  for 
exai^le,  the  long  run  proportion  of  time  one  of  the  servers 
is  idle,  the  average  expected  busy  period  of  one  sesvor*  and 
the  average  tioM  it  takes  a oustoaer  to  complete  one  cycle 
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of  the  network.  There  is  some  indication  (cf.  Lewis  and  Shedler 
(1973])  that  the  service  times  for  one  of  the  servers  should 
be  correlated.  One  object  of  this  paper  is  to  show  that 
correlation  of  the  service  times  does  make  a difference  in  the 
limiting  behavior  of  the  network. 

Let  S (respectively  X_)  denote  the  nth  service 
time  for  server  2 (respectively  server  1).  There  are  many  ways 
in  tdiich  to  use  EARMA(1,1)  processes  to  model  dependence  within 
each  sequence  of  service  times  and  cross  correlation  between 
the  two  sequences.  In  this  paper  we  will  assume  that  {S^} 
is  a sequence  of  independent  random  variables  each  with  an 
exponential  distribution  with  mean  0 < y < <»  and 

is  am  EARMA(1,1)  process  independent  of  iS^}  with  mean 
0 < X < «,  and  parameters  6 and  p.  Other  queueing  models 
using  EARMA  processes  will  be  considered  elsewhere. 

Let  Z(t)  denote  the  number  of  customers  both  waiting 
and  being  served  by  server  1 at  time  t.  The  process 
Z » {Z(t);  t ^ 0}  takes  the  values  {0,1,..., N}.  Let  W(t) 
be  the  virtual  waiting  time  at  server  1 at  tisw  t;  that  is, 
lf(t)  is  the  sum  of  the  service  times  of  the  customers  in  queue 
at  server  1 at  time  t and  the  remaining  service  time  of  the 
customer  currently  being  served. 

In  the  next  two  sections  %»e  will  obtain  limiting  results 
for  l(t>  and  1f(t)  aa  t • for  the  case  N ■ 2.  wa  will 
idaoir  that  the  limiting  distribution  of  I (t)  in  the  BARNA 
ease  is  the  sams  as  in  the  ease  in  whicdi  the  service  times  of 
sarmar  1 arm  indepwidaoti  this  result  is  also  true  for  the  cms 
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N > 1 and  seems  to  be  due  to  the  constraints  Imposed  on  the 
net%fork  by  having  only  one  or  two  customers.  The  long  run 
average  virtual  waiting  time  is  then  computed  for  the  case 
N * 2 and  is  found  to  be  different  from  the  independent  case. 
This  is  because  the  positive  correlation  between  the  service 
times  of  server  1 tends  to  Increase  the  virtual  waiting  time. 

In  Section  4 we  analyse  the  network  for  N > 2 customers 
and  show  the  existence  of  a limiting  distribution  for  Z(t) 
as  t •*■  ».  We  then  present  some  simulation  results  for  the 
limiting  distribution  of  Z(t)  as  t -«■  <»  for  the  case  N « 5. 
The  results  show  that  the  correlated  service  times  do  make  a 
considerable  difference  in  the  limiting  distribution. 


2.  THE  TWO  CUSTOMER  CYCLIC  QUEUE 

In  this  section  we  will  stt^y  the  process  Z « {Z(t);t  ^ 0} 
end  obtain  the  limiting  distribution  of  Z(t)  as  t for 

the  queueing  network  of  Section  1 with  N ■ 2. 

2.1  Prelisdnaries 

We  define  here  three  seqiMnces  %diioh  will  be  needed 
in  the  analysis.  For  concreteness  we  will  elways  arausM  that 
Z(0)  « 0.  Let  be  the  inareasing  sequence  of  arrival  and 

dsparture  tiaws  for  server  1 that  includes  all  departure  tisMs 
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and  those  arrival  times  at  server  1 that  occur  when  server  1 
is  idle.  Let  be  the  number  of  customers  both  waiting 

and  being  served  just  after  time  T^.  More  precisely  put 
Yq  = Z(0)  =0  and  = 0 and  define  Y^,  n ^ 1, 
recursively  as  follows.  If 

Vl  - = Vl  - -1- 

If  Y_  > 0,  let 
n 

(2.2)  = inf{t  > T^:Z(t)  < Z(t-)}  , = Z(T^^j^+0)  . 

For  exaoqple,  since  Z(0)  * 0,  is  the  time  of  the  first 
arrival  of  a customer  to  server  1 and  Y^  is  the  number  of 
customers  at  server  1 just  after  time  which  must  be  1; 

is  the  time  of  completion  of  the  first  service  for  server  1 
and  Y2  is  the  number  of  customers  both  waiting  and  being 
served  at  server  1 just  after  time  T^;  the  number  could  be 
G or  1.  In  general,  if  ^ then  the  time  of 

ccmipletion  of  the  next  service  for  server  1 after  time  T 

n 

and  the  nundder  of  customers  both  waiting  and  being 

served  at  server  1 just  after  time  If  * 0,  then 

Vl  is  the  time  of  the  next  arrival  to  server  1 and  Y^^^  « 1. 
Mote  that  since  we  are  consitering  the  case  in  which  M « 2, 

Y^  can  <Mily  take  the  values  0 and  1. 

Let  denote  the  autoregressive  part  of  the  next 
senrioe  to  he  eoaqpleted  after  time  More  precisely,  put 

Lg  • 0 and  recursively  define  * inf{n  > h|(<Y^  > 0); 
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0 . Then  since 


L.  is  the  index  of  the  kth  time  that  Y “ 

K n 

Z(0)  * 0,  Aq  * Aq,  = Aq  and  A^  = ^n-k-1 
L(k)  < n _<  L(k+1)  . 

Since  the  service  times  at  server  1 are  dependent  random 

variables,  the  dependence  of  {2(t);  t > T } on  {Z(t);  t < T } 

n — n 

is  not  only  through  but  also  through  the  service  times  that 

were  completed  before  time  However,  by  (1.1)  and  (1.2), 

{Z(t);  t > T^}  is  conditionally  independent  of  {Z(t);  t £ T^} 

given  (A„,Y„) . Further  {(A^,Y^)}  is  a discrete  time  Markov 
n n n n 

process  with  state  space  (»_^  x {0,1},  a (g_|^  x {0,1}))  where 
* (0,“>],  denotes  the  Borel  subsets  of  and 

0(8^  X {0,1})  denotes  the  product  a-algebra  generated  by 
and  the  subsets  of  {0,1}. 

Let 

^j(*'B)  - « B,  - J|A„  - X,  - i) 

and  P(x,B)  be  the  matrix  whose  i,j  entry  is  P^^(x,B) 

for  i,  j € {0,1},  X € and  B € g^.  The  transition 

probabilities,  P^j(x,B)  are  easily  derived.  If  Y^  ■ 0, 

then  tine  of  the  next  arrival  after  time  T^ 

and  hence  Y^^j  • 1 and  * *n*  ^n  “ ^ there  Is 

an  arrival  during  the  service  tine  Vl  - '^n'  then  Vl  - 

If  there  Is  no  arrival  during  the  service  tine  T_.,  - T , 

n+i  n 

then  ■>  0.  Since  {X^}  Is  an  EARNA(1,1)  process  and 

{8^}  Is  a sequence  of  Independent  exponential  random  variables 
with  sMsn  %»e  have  the  following  swtrlx 
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0 


(2.3)  P(x,B) 


6 (B) 
X 


-wx 

Ete  ^;Aj^€B|Ajj=xl 


EU-e 


-ViX, 


;Aj^€B1Aq=x] 


where  is  1 if  x € B and  0 otherwise. 

We  will  now  define  other  processes  related  to  {A  ,Y  ,T  } 

n n n 

which  will  also  be  used  in  the  analysis.  For  T < t < T 

■*  n — n+1' 

let  Y(t)  = Y^,  A(t)  = A^,  and  U(t)  = t-T^.  Note  that  Z{t)=0 

if  and  only  if  Y(t)  = 0;  Z(t)  =1  if  and  only  if  Y(t)  =1 

and  there  is  no  arrival  at  server  1 in  the  time  interval 

(t-U(t),tJ;  and  Z(t)  -2  if  and  only  if  Y(t)  =1  and  there 

is  an  arrival  in  the  interval  (t-U(t),t].  Hence  the  limiting 
behavior  of  Z(t)  as  t ^ is  related  to  that  of 
(A(t) ,Y(t) ,U(t) ) . Furthermore,  the  limiting  behavior  of 
(A(t) , Y(t) ,U(t) ) depends  on  that  of  (A^,Y^) . As  a result, 
we  will  first  compute  the  limiting  distribution  of  (A^,Y^); 
then  use  it  to  confute  the  limiting  distribution  of 
(A(t) ,Y(t) ,U(t) ) ; and  finally  compute  the  limiting  distribution 
of  Z(t)  as  t 

2.2  Limiting  Properties  of  ’ 

Fix  a Borel  subset  B of  K ^ with  positive  Lebesgue 
measure.  Prom  (1.2)  it  follows  that  (A^)  is  a discrete  time 
Harkov  process.  To  show  that  the  process  is  recurrent  note  that 
expression  (6.2)  of  Jacobs  and  Lewis  [1977]  for  the  kth  order 
transition  probability  of  {A^>  ingtlles  that  for  each 
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b > 0 and  0 < 6 < / Xe~  ^ dx  there  exists  a k such  that 
“ B 


Hence 


inf  Q (x,B)  > « 
x€  [0,bl 


P(  U € B}  IAq  - X)  > 6 , 

ni*l 


X € K, 


Therefore,  by  Proposition  (5.1)  of  Orey  [1971] 


P(Aj^  € B infinitely  often  |Aq  * x)  ■ 1 

for  all  X and  thus  {A^}  is  recurrent  with  respect  to 
Lebesgue  measure  in  the  sense  of  Orey  [1971,  page  4]. 

Since  the  service  times  for  server  2 are  independent 
with  comoon  exponential  distribution,  a similar  argument  shows 
that 


(2.4)  ^ ^n  “ ^ infinitely  often  | A^-x,  YQ“i)  ■ 1 

for  i»  j € {0,1},  X € . Thus,  by  Theorem  (7.1)  of  Orey 

[1971]  there  exists  a possibly  o-fiuite  invariant  measure 
{Vj(B);  j * 0,1,  B € for  the  transition  probability  of 

((A„,Y„)). 

The  invariant  measure  satisfies  the  following  system 


of  egtshtions 


/ ir-Cax)  Pq  q(x,B)  + / ir.(dx)  P.  o(x,B)  - 

qv  «»v  q*  u 


Xg(dx)  Pq^j^(x,B)  + / ffj^Cdx)  “ w^(B)  . 


Rewriting  the  two  equations  using  the  matrix  (2.3)  yields 


-MX, 


(2.5)  / ir^(dx)E[e  "‘■j  € B|Aq»x]  « ^q(B) 


« -MX, 

(2.6)  Xq(B)  + / irj^(dx)  EU  - e ■‘•jAj^  € B|Aq»x1  » iTj^(B)  . 

Stibstituting  the  expression  for  ^^^(B)  from  (2.5) 
into  (2.d),  equation  (2.6)  becomes 

00 

(2.7)  / Trj^(dx)  P(Aj^  € B1Aq»x)  « . 

By  the  result  after  (6.2)  in  Jacobs  and  Lewis  [1977],  equation 

(2.7)  implies  that 

(2.8)  W| (B)  ■ c / Ae“^^  dy 

* B 

for  any  ncm-negative  constant  c.  Substituting  the  expression 
for  V,  into  (2.5)  we  have 


(2.9) 


« **)jX 

iTq(B)  ■ C E(e  € B|Aq-x]  dx  . 

He  now  want  to  choose  c so  that  ir  is  a probability  and 
hence  ffj|^(B)  » € B,  To  this  end  we  set 

1 » 1^1  (»+  ) + ^ 

- c(l  + / Xe"^*  B(e  ^|Aft-xJ  dxj 

0 “ 

- c(l  + X(X+w)“^] 

since,  if  A^  has  an  exponential  distribution  with  mean  X~^, 
so  does  X^.  It  now  follows  that  c * (X-fy)  (2X<fy) 

This  result  will  be  used  in  the  next  subsection  to  compute  the 
limiting  distribution  of  (A(t) ,Y(t} ,U(t) ) . 


2. 3 Limiting  Properties  of  an  Imbedded  Semi -Markov  Process 

Since  {(A„,Y.),T„}  is  a Harkov  renewal  process  in  the 
n n n 

sense  of  ginlar  (1975],  { (A(t)  ,Y(t)  ,0(t) ) ;t  ^0}  is  a Markov 
process.  Hence,  { (A(t)  ,Y(t) ) it  ^0)  is  a semi-Markov  process 
of  the  second  type  in  the  sense  of  Jacod  (19731 . He  will  use 
the  results  of  Jacod  (19731  to  ooaipote  the  limiting  distribution 
of  (Aft)  ,T(t)  ,u(t) ) as  t iMiioh  will  then  be  used  to 
compute  the  limiting  distribution  of  k(t). 

Since  (e^l  and  {8^}  are  indapendant  saquanoaa  of 
indapandant  axponantial  variables,  the  prooaas  {A(t),y(t)} 
is  right  ccHitinuoua  with  left  hand  liadta  aatiafying  tha 
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hypothesis  R-3  on  page  85  of  Jacod  [1973 ] concerning  the  set 
of  discontinuity  points  of  (A(t),Y{t)).  Furthermore, 

{ (A(t) ,Y(t) ) ;t  ^ 0}  is  recurrent  in  the  sense  that,  if 
B € has  positive  Lebesgue  measure,  then 

00 

/l^jj«Y(s)  lg*A(e)ds  - -|Aq  - X,  Yq  - i}  - 1 

for  i,  j i {0,1}  and  x € where  lg(x)  is  1 if  x € B 

and  0 otherwise.  Therefore,  Theorem  III-IO  on  page  103  of 
Jacod  [1973]  applies  to  show  that 

00  00 

(2.10)  Vq(BxC)  » d / ir-(dx)  / e"^'®  l_(x)  l„(s)ds 

V 0 0 B c 

«»  00 

(2.11)  v^(BxC)  - d / irj(dx)  / P{X^>  s|Aq«x}  lg(x)  l^(s)ds, 

B,  C € 2^,  is  an  invariant  measure  for  the  transition  function 
of  the  Markov  process  {(A(t),Y(t) ,U(t)) ;t  > 0)  for  any 
d > 0. 

He  will  now  choose  d so  that  v^,  i - 0,1,  is  a prob- 
ability. Put 

1 - v^(*J)  ♦ Vj(»J)  - dlii*^ 

- de|X[|i(A4>M)l‘^  ♦ . 


BoBea. 


■Pi^  (2.12)  de  - Xw(X+y)  (X^+XiH-y^)"^ 

and  d ■ Xy(X+2y) (X^+Xy+y^) 

2.4  The  Limiting  Distribution  of  2(t)  ae  t * 

Using  (2.10),  (2.11),  and  (2.12)  we  cem  now  compute  the 
limiting  distribution  of  Z(t),  the  nusiber  of  customers  both 
waiting  and  being  served  at  queue  1,  as  t -»■  <».  Recall 
that  we  are  considering  the  case  N * 2 and,  therefore,  Z(t) 
has  the  values  0,  1,  or  2. 

From  the  argument  after  (2.3)  it  follows  that 


(2.13)  lim  P{Z(t)-0|Z(0)-0}  - lim  P{lf(t)-0 (y (0)-0}  » Vq(b;) 
t-*--  t ^ 

» X^(X^  + Xy  + y^l“^; 


i : 


(2.14)  lim  P{z(t)-l|Z(0)-0)  - iim  rt(t>*i|y(oy-dJ 

t -►  OB  t • 


I f v^fdx^ds)  e 


' Xyix^  + Xy  + y^r^ 

; } " ( " Hi  ^ ^ ■*  7 


by  (2.11)  and  (2.12);  and  finally 


f.  + ' t ( - ■»  i 'O  ' ■ 
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II 


(2.15)  li«  P{a(t)-2|Z(0)-0}  - li»  Y<0)-0J 

t *►  « t • 

« y^IA^  + Ay  y^j"^  . 

Bote  that  the  Uniting  distribution  £or  the  nunber  of 
custoners  waiting  or  being  served  by  server  1 with  EARMA(1,1) 
service  tines  is  the  sane  as  if  the  service  tines  for  server  1 
were  independent  randon  variables  each  with  an  exponential 
distribution  with  nean  a"^  (cf.  Gaver  and  Thonpson  [1973]). 
ffe  feel  that  this  is  due  to  the  constraints  inposed  on  the  net- 
work by  having  only  two  custoners.  This  conjecture  is  shown 
to  be  true  by  the  simulation  results  of  Section  4 %fhere  it 
beccmes  clear  that  the  result  is  not  true  if  N > 2* 


3.  TOE  VIRTUAL  WAITING  TI»e  FOR  THE  CYCLIC  QUEUE  WITH  TOO 

CUSTOMERS 

In  this  section  we  will  compute  the  long  run  virtual 
waiting  tine  for  server  1 in  the  case  in  %rhich  the  service 
tines  for  server  1 form  an  BARNA(1,1)  process  and  there  are  two 
ctistoners  in  the  systen. 

The  virtual  waiting  tine  for  server  1 at  tisw  t,  W(t)* 
is  the  sun  of  the  service  tines  of  the  custosMrs  in  queue  at 
server  1 at  tiaw  ,t  plus  the  regaining  service  tine  of  thm 
custoswir  currmtly  being  served,  to  define  N(t>  nore  precisely. 


i 


i 


I 

let  N(t)  - 8up{n:T^  5.  t}.  1£  Z(t)  la  1 or  2,  then 

V(t)  ■ ^ timm  from  t until  the  completion 

of  the  current  service.  If  Z(t)  - 2,  then 

is  the  length  of  the  service  tlsie  for  the  custoner  who  arrived 
at  server  1 in  the  time  Interval  (t-U(t),t].  Hence,  the 
virtual  waiting  time  at  time  t Is  defined  formally  as 


(3.1) 


W(t) 


I ° 

I V(t) 

y(t)  + “ ’^N(t)+l 


if  Z(t)  - 0 , 
if  Z(t)  - 1, 
If  Z(t)  - 2 . 


Prom  (3.1)  it  follows  that  { (A(t) ,y (t) ,Z(t) ,U(t) ,V(t) ,W(t)) » 
t ^ 0}  Is  a Mar)cov  process.  Further,  the  process  Is  recurrent 
in  the  sense  of  Azema,  Duflo,  and  Revuz  [1969].  Thus  by  the 
ergodic  theorem  in  Section  (3.1)  of  that  paper,  there  exists  a 
constant  W such  that 

^ 1 t 

(3.2)  W - lim  i / W(s)d8 

t ^ 0 

P(*|A(0)«x,  Z(0)«0) -almost  surelif  and 

(3.3>  9f  • lim  ^ 11/^  1»(a)ds|  A(0)*X,  Z(b)-O) 

t • 0 

for  almost  ell  x.  Me  will  oee  (3.3)  to  compote  9. 

Pixet,  fpom  (3.11  and  the  axgoileat  after  (2.3) 
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I 

\ 

i 

I 

i 

I 


I 


(3.4)  B(1f(s)  |A(0)-x,  Z(0)-0] 

- » y(8)-l|A(0)«x,  2(0)»0] 

- EIV(b);Y(*)-1|A(0)-x;  Z(0)-01 

tfa  will  first  coapote 

(3.5)  li«  i / E(V(8)j  Y(8)*1|a(0)  - X,  Y(0)  - 01  ds  . 

t ' 0 

Sines  (A(t)  ,Y(t)  ,U(t) ) has  a liaiting  distribution  as  t 
by  ths  proof  of  Theorem  III-IO  on  page  105  of  Jacod  [1973], 
we  have  that 

U, (D)  - 11m  P{Y(t)-l,  (A(t) ,0(t) ,V(t) ) € DlA(0)«x,Y(0)-0} 

^ t ♦ • 

exists  and  equals 

**  -1  **  y 

cd  f Xe"**  ds  / P{X,  € dy)Ag-s)  / dx  l_(s,x,y>x) 

0 0 0 

where  D is  a Borel  subset  of  . By  the  ergodlo  theorem  of 
As^ma,  Dttflo  and  Rsvus,  (3.5)  fmuals 

(3.5)  / |i^(ds,dy,ds)s  • ii(X«»-|i)  Cl(X*+Xy+»i*)  1*^  . 

By  (1.1)  and  (1.2) » the  second  tens  of  (3.4) 


(3.7)  E(  *(0)-0] 

y(8)-l|A(0)-x,  Z(0)-0J 

+ (l-e)E(pA(8)(l-«“^^^*^);  y(«)-ljA(0)«x,  Z{0)-0J 
+ (1-8)  (1-p)  Ete^(gj^^(l-e"*'^^®^;Y(s)-l|A(0)-x,Z(0)-01. 
Since  e))(g)4.2  Independent  of  (y(e),U(s)) 

(3.8)  lim  i 

• A ^ lia  f f P{Z(s)*2iA(0)»x,Z(0)»0)ds 

t ^ 0 

- V^(A(A^+Au+u^)]"^ 

by  (2.15)  end  the  previously  cited  ergodic  theorem. 

The  ergodic  theorem  elso  implies  that 

(3.9)  lim  y B(A(s)  (l-e’*'”^*^);y(8)-l|A(0)-x,Z(0)-01d8 

t * • ^ h 

■ / V,  (ds,dx)t(l-e"*'*) 

- cd  / Ae*^S  «s  ^ P(X^  > x|A^-t}  (l-e“»*“)dx 


cd  / Xe"***  dz  / ieP{Bej  > x}+(l-6)p{Bej^>  x-z)]  (l-«’’*'*)dx 


(X{X^+Xu+u^] (X+u) (X+Bm)!"^ 


xy^t3X^+2Xu+B(2Xy+2y^-X^)  + B^ (-2Xy-y^-X^) J 


by  (2.11),  (1.1)  and  (1.2). 

tie  now  ttirn  our  attention  to  the  confutation  of 

(3.10)  ^lim  i ECejj^^j^^(l-e"‘'"^®^,y(8)-l|A(0)-x,2(0)-0]d8  , 

First,  note  that  ^^n'^n'^ntl^  * discrete  tine  Marlcov  process 

with  limiting  distribution  ir^(Bxc)  ■ Zj^(B)  /j,  Xe'^^dx  for 
B,  C € R^.  Hence,  by  similar  arguments  to  those  in  subsection  2.3, 
(A(t)  ,Y(t)  ,e|j^^^^^,U(t) ) has  a limiting  distribution 
Further,  for  B,  C,  D € g^,  by  Theorem  III-IO  of  Jacod 

"v.  (BxCxo)  ■ d/  z.  (dz)  /Xe”^*dx  / P(x  > s|A„*z,  e.Bx)ds 
* B*C  dI''-'- 

- d /it,  (dz)  / Xe’^*dx  /(Bl|.  (Bx)-»-(1-B)1-  (Bx+z)}d8 

B C D ^ 


le  constant  d by  (1.1)  and  (1.2). 
To  evaluate  note  that 


v^(Bxd)  - ^j(1  * x D) 


I 


H«nctt,  from  (2.8)  and  (2.11),  d > d.  By  the  ergodlc  theorem 
again,  (3.10)  equals 


I 

\ 

i 


(3.11)  / V,  (dz,dx,ds)  x(l-e*'^*) 

- (X(X^+Xy+p^) (X+U6)]“^  y^(X(l+B+$2)  + y(B+e*)l 

after  some  slnqpllflcation. 

Putting  together  (3.3)-(3.11)  we  obtain  after  some 
slnq>llfication 

(3.12)  IX(X^+Xy+y^) W 

« yX  ♦ 2y^ 

+(1-8) y^ [ (X+y) (X+By) ]"^{ (8X(X+y)+8^ (X+y)^) 

+p I2X (X+y) +B{-2X*+y2) .202 ( x+y) } . 

Putting  Y * the  traffic  intensity,  (3.12)  becoows 

(3.13)  X?i  - (2+Y)(Y^n+l)“^  + (l-B)C(Y*+Y'*'l)(Y+l>(Y+a)l*^ 

» { I By  ( Y-M  ) ♦B^  ( Y+1)  *1  +0  f2Y  (Y+1 ) ♦B  (-2y^+1  ) - 2B*  ( Y+D  ^ 1 ) • 

...  - , . _ V 

If  B • 1»  then  til#  serwiee  tUmB  for  servor  1 are 
independent  exponsntiel  nendofi:  wes^dWLelr  Ifithi  Mini  and 

rf  y.; 

C2.14)  xii^4a4  - (2+tHY*4?|(|i-l)*^ 


i 


as  expacted  in  this  case.  Note  that  W ^ and,  hence, 

the  positive  correlation  of  the  service  tines  increases  the 
average  virtual  waiting  time.  Also  there  is  a 0q  such  that 
for  0 < Increasing  p increases  while  fcr  0 > 0q 

increasing  p decreases  W.  The  value  0q  . is  a solution  to 
the  quadratic  equation 

(3.15)  -2(y+1)^0^  + (-2y+1)0  + 2y(y+1)  « 0 

and  hence  depends  on  the  traffic  intensity  Y . 

To  give  an  idea  of  the  effect  of  different  0 and  p 
on  W we  give  Table  1 whose  entries  are  differences  between  ^ 

and  £or  various  values  of  0 and  p for  » 1 

in  this  case) . 


p\0 

0 

.1 

.3 

.5 

.7 

.9 

0 

0 

.03 

.09 

.11 

.10 

.04 

.1 

.05 

.07 

.10 

.11 

.09 

.04 

.3 

.15 

.15 

.14 

.12 

.08 

.03 

.5 

.25 

.23 

.18 

.13 

.08 

.03 

.7 

.35 

.30 

.21 

.13 

.07 

.02 

.9 

.45 

.38 

.25 

.14 

.06 

.01 

TlltLI  1. 

Values 

of  « 

•i-d 

for  the 

: 

m 

o 

X • p ■ 1 

Note  that  it  is  the  autoregressive  part  of  the  EARM2l(l,l) 


service  time  that  causes  the  most  change  in  W.  As  expected, 
for  3 small  W increases  with  increasing  p;  for  0 large 
W decreases  with  increasing  p;  the  value  at  which  the  change 
occurs,  3q,  is  .55  in  this  case. 


4.  THE  CYCLIC  QUEUE  WITH  A FINITE  NUMBER  OF  CUSTOMERS 

In  this  section  we  will  study  the  queueing  network  of 
Section  1 in  the  case  in  which  there  are  N ^ 3 customers.  We 
will  use  without  mention  previous  notation  adapted  to  the 
present  case.  As  before,  we  are  interested  in  the  limiting 
distribution  as  t ^ of  Z(t),  the  number  of  customers  either 
waiting  or  being  served  at  server  1 at  time  t.  Again,  we  first  ■< 

consider  the  discrete  time  Markov  process  which  now 

has  state  space  ( ^ tS,l,*.. ,N*1 J , 0 ( ^ £0,1,..., N*1 J ) ) 

where  9 x {0,1, . . . ,N-1})  denotes  the  product  a-algebra 
of  the  Borel  subsets  of  and  the  subsets  of  {0,1, . . . ,M-1} . 

L*t  P^j{x,B)  - € B,  »„+i-j|A„-x,  Yjj-i)  for  B € g^, 

X € , and  i,  j € {0,1,.. If  P(x,B}  denotes  the 

matrix  whose  (i,j) -entry  is  P^j(x,B),  then  it  is  not  hard  to 
show  that 


" 0 

5 (B) 

0 . . 

. . . 0 

0 

bQ(x.B) 

X 

b^(x,B) 

bj  (x,B) 

^N-2 

aj^l(x,B) 

0 

bQ(x.B) 

bj^(x,B) 

bj,_3(x,B) 

0 

0 

bQ(x.B) 

bjj_^(x,B) 

ajj_3(x,B) 

P(X,B)  - 

• 

• 

• 

• 

• 

• 

• 

• 

• • 

• 

• 

• 

• 

• • 

• 

0 

0 

0 . . 

bo(x,B) 

• • • u 

aj^  (x,B) 

where 

r -uX  (MX  I 

b,^(x,B)  -Be  ^ € BjAp-x 

and 

m 

a.  (x,B)  = I b (x,B)  . 
n-k  " 

By  similar  methods  to  those  of  subsection  2.2  one  can  show 
that  there  exists  a possibly  a-£inite  invariant  measure  for  P 
which  satisfies  the  following  system  of  equations  for  B € 


(4.1)  / ir^(dx)bQ(x,B)-irQ(B) 


/ Wp(dx)4jj(B)-*'/  ir^(dx)b^(x,B)<«'/  X2(dx)bQ(x,B)«v^(B) 


(4.2) 


(4.3)  / irj(dx)b2(x,B)+  / ^2  (dx)b^  (x,B)+  / tdx)bg  (x,B)-X2 (B) 


(4.4) 


/ Wj^(dx)ajj_j^(x,B)+/ff2  (dx)ajj_,2t^»®^‘*’* 


’^N-l 


(dx)a^  (x,B)»Xjj_ 


Substituting  the  expression  for  ifQ(B)  of  (4.1)  into  (4.2) 
and  then  adding  all  the  equations  except  (4.1)  together  yields 
the  equation 
00 

^ (dx)  P{Aj^€  B|A^»x}  ■ Iirj^+ir2+*  (B) 

By  the  result  after  (6.2)  of  Jacobs  and  Lewis  [1977]  we  have 

+•••+  - c»  / Xe*^*  dx 

B , 

for  any  constant  c*.  Hence/  we  can  and  will  choose  c*  so 
that  the  invariant  measure  x is  a probability* 

We  will  now  consider  the  limiting  distribution  of  2(t) 
as  t . Let  K(t}  be  the  nuaber  of  service  ccxapletions  by 

. ' « ■ ! , i . i;. ' ' ‘ - i . -x  ! , ■ ‘ '1  . 

server  2 in  the  time  interval  it  - U(t)  /t] . Binoe  the  seirvioe 
times  of  server  2 are  independent  each  with  the  same  enponmtial 
distribution  {(A(t)/y(t)/0(t}/X(t})it  ^0}  is  a Narlcov  prooess. 
Note  that  {t(t)»0}  » {y(t)«0}| 


(B) . 


22 


and 


1 

{Z(t)-i}  - U {Y(t)-k,K(t)-i-k}  , 
k«l 

N-1 

{Z(t)«N}  » U {Y(t)-k,K(t)-N-k}  . 
k-1 


0 < 1 < N; 


Thu8»  using  the  techniques  of  subsections  2.3  and  2.4  one  can 
show  that 


1-1. 


v<0)  - Urn  P{Z(t)-0|Z(0)-0,A(0)-0}  • dy  t 

t 


v(H)  « lin  P{Z(t)-N|Z(0)«0,A(0)-0} 
<0 

H-1  « 


^ I I if,.(dy)  I P{X->s|A-y}  I e"^*  ^ dsy 
k-1  6 0 * * j-H-k 


j 


and 


v(i)  - lin  P{Z(t)-iiZ(0)-0,A(0)-x} 

t > • 

i • 


■ ^ Jl  9 0 '‘*1**I*0*J'’ 


i-k 


ds 


for  soipe  positive  constant  2. 

It  ssews  difficult  to  solve  (4.1)*(4.4}  for 
i - Of...,ll-l.  Bence,  «fe  are  unable  to  obtain  excpliolt  expressions 
for  the  liadting  distribution  of  t(t)  as  t as  we  could  in 
the  ease  B»2.  lie  will,  however,  give  eoas  siaiiIat4.on  results 
to  ia^oate  the  lijsiting  behavior  of  S<t)  as  t for  hif^imr 
veiess  of  «. 

m siasaation  is  based  oe  tlMi  fellewinf  observa^on^ 

•iebs  i»e  servioe  tiass  of  eervsr  t ete  iadaystideet  and  sspeeewtthlly 

tJ 


distributed  { (A(t) ,Y(t) ,U<t) ,Z (t) ) ;t  ^0}  is  a Markov  process 
that  satisfies  the  hypotheses  of  the  ergodic  theorem  of  Azfoa, 
Duflo  and  Revuz  [1969].  Hence, 


1 ^ 

lim  -z  I l|.i«Z(a)ds  « v(i)  , 
t ^ 0 


1 - 0,...,N 


almost  surely  P(>|A(0)«x,  Z(0)«0).  Further,  if  is  the 

time  of  cong>letlon  of  the  Mth  service  time  for  server  1,  then, 
since  « as  M -*-  « almost  surely  with  respect  to 

P(*  |A(0)-x,Z(0)-0) 


(5.5) 


1 " 

lim  — / lr,i*Z(s)ds  ■ vd),  i ■ 0,...,N 

M 0 


almost  surely  with  respect  to  the  same  probability. 

The  following  results  are  from  a simulation  designed 
by  Professor  P.  A.  H.  Lewis  at  the  Naval  Postgraduate  School, 
Monterey,  California  for  a closed  queueing  network  having 
BARMA(1,1)  service  times,  nie  simulation  we  used  consists  of 
the  computation  of 


(5.6) 


1 /m 


for  If  equal  to  the  SOOOth  service  time.  Thm  computation  was 
repeathtf  for  100  iiwwpendent  replieatioos  and  the  sample  mean 

- - Si,  -Vii  .t-vS  ^ ^ I 

and  iriurisnee  over  the  lOO  replieations  were  oompiKted.  All 
rittis  were  performed  cm  an  ZM  360/67  oommitor  at  tte 

Naval  fi^iigra^ta  lo^l  uai^  the  Zd4HlilOQN  padkife  (fcaawneiiih  . 


and  Lewis  [1973] ) which  generates  numbers  according  to  the 


schesw  given  by  Lewis,  Goodman,  and  Miller  [1969]  and 
exponentially  distributed  random  numbers  using  the  Marsaglia 
"rectangle-%fedge-tail"  method.  Tests  of  the  random  number 
generator  are  given  in  Learmonth  and  Lewis  [1974]. 

Table  2 gives  the  results  of  a simulation  of  the  queueing 
network  for  the  case  H«5  and  X«ii*l.  The  entries  in  the  t2U>le 
are  the  sample  means  of  (5.6)  over  the  independent  replications 
for  different  values  of  6 and  p.  The  nusdbers  in  the 
parentheses  are  estimates  of  the  standard  deviations  of  the 
estimate  of  v(i).  For  comparison,  note  that  the  limiting 
distribution  in  the  case  in  which  server  1 has  independent 
exponential  service  times  and  X«y*l  is 

(5.7)  lim  P{Z(t)-i)  « 1/6  - 0.1666,  i - 0,...,5  , 

t * 

(cf.  Gaver  and  Th<mq>son  [1973]).  This  corresponds  to  the  case 
6 > 0,  p « 0 in  the  first  line  of  the  table. 

Note  that  again  it  is  the  autogressive  part  of  the 
service  times  that  causes  most  of  the  change  in  the  estimates 
for  the  limiting  distribution  of  Z(t)  as  t . The 
positively  correlated  service  tisms  increase  the  probabilities 
of  server  1 or  server  2 being  idle;  they  also  increase  the 
probability  of  all  cestoasrs  being  in  one  or  the  other  service 
center.  This  seems  to  be  due  to  the  fact  that,  if  p is  large 
and  8 is  small,  tbmn  having  a large  servios  time  at  on*  time 

2S 


1 

0 

1 

2 

3 

4 

5 

8-  0 

pm  0 

.17(.009) 

.17(.007) 

.17(.005) 

.17(.006) 

.17(.007) 

.17(.009) 

.1 

.18(.01) 

.16(.007) 

.16(.005) 

.15(.00S) 

.17(.007) 

.17(.01) 

.5 

.21(.01) 

.13(.006) 

.12(.006) 

.12(.005) 

.18(.007) 

.21(.01) 

.9 

.26(.02) 

.ll(.Ol) 

.079 (.006) 

.088(.006) 

.20(.01) 

.26(.03) 

B-.l 

pm  0 

.17(.01) 

.16(.008) 

.16 (.005) 

.16(.005) 

.17(.007) 

.17  (.01) 

.1 

.18(.01) 

.16(.008) 

.16(.005) 

.15(.005) 

.17(.007) 

.17(.01) 

.5 

.21(.01) 

.15(.007) 

.13(.006) 

.13(.005) 

.18(.007) 

.20(.01) 

.9 

.25(.02) 

.13(.01) 

.092 (.006) 

.10(.005) 

.18(.01) 

.25(.03) 

6-. 5 

p-  0 

.18(.01) 

.17(.007) 

.16 (.005) 

.15(.006) 

.16 (.007) 

.18(.01) 

.1 

.18(.01) 

.17(.007) 

.16 (.005) 

.15(.006) 

.16(.007) 

.18(.01) 

.5 

.19(.01) 

.17(.007) 

.15(.006) 

.15(.006) 

.16(.007) 

.18 (.01) 

.9 

.20(.02) 

.16 (.01) 

.14(.007) 

.14(.006) 

.16(.01) 

.20(.02) 

8-.  7 

pm  0 

.18(.01) 

.17(.007) 

.16(.005) 

.15 (.005) 

.16 (.007) 

.18(.01) 

.1 

.18(.01) 

.17(.007) 

.16(.005) 

.15(.006) 

.16(.007) 

.18(.01) 

.5 

.I8(.01) 

.17(.008) 

.16 (.005) 

.15 (.006) 

.16(.008) 

.18(.01) 

.9 

.18(.01) 

.17(.009) 

.16(.006) 

.16(.006) 

.16(.009) 

.17 (.01) 

6-.9 

p-  0 

.17(.01) 

.17(.007) 

.16 (.005) 

.16 (.006) 

.16(.007) 

.17(.01) 

.1 

.17(.01) 

.17(.007) 

.16(.005) 

.16(.007) 

.16(.007) 

.17(.01) 

.5 

.17(.01) 

.17(.007) 

.16(.005) 

.16(.006) 

.16(.008) 

.17(.01) 

.9 

.17(.0l) 

.17(.008) 

.17(.005) 

.16(.006) 

.16(.008) 

.17(.01) 

TABLE  2.  Estiaaccs  for  th«  Halting  distribution  for  tbs  tambrnr  of  custoasrs 
In  qusue  for  N"S  custoasrs  and  ssrvur  1 has  EAIDfA(l,l) 

ssrvics  tiass  and  aaruar  2 lias  iada|Madaot  aapoaantlal  sarvlca  elasa. 


•fL- * ll 


implies  that  the  service  times  will  tend  to  be  large  for  a 
while;  having  a small  service  time  implies  that  the  service 
times  will  tend  to  be  small  for  a while.  The  dependent  service 
times  appear  to  have  little  effect  on  the  long  run  average 
number  of  customers  that  are  waiting  or  being  served  at  server  1; 
the  average  is  2.5  for  all  values  of  (3/P)  except  the  values 
(0,.9),  (.1,.5),  and  (.1,.9)  for  which  it  is  2.6.  This 
invariance  reflects  the  fact  that  the  chimge  in  the  limiting 
distribution  {v(i)}  due  to  the  dependence  is  somewhat 
symmetric  in  i . 

Finally,  %«e  present  in  Table  3 the  results  of  a simulation 
to  investigate  what  happens  if  the  service  times  of  both  servers 
are  independent  EASMA(1,1)  processes.  In  that  simulation 
tne  nundser  of  customers  N is  5 and  the  parameters  for  the 
EARMh(l,l)  service  times  for  server  1 are  3^  » .1  and  ■ .9. 
The  parameters  ^2  BARNh(l,l)  service  times 

for  server  2 are  allo%«ed  to  take  on  several  values.  The 
expression  (5.6)  is  conqputed  for  M > 10,000  for  750  independent 
realisations.  The  entries  in  the  table  are  the  values  for  this 
sisnilation  of  the  same  quantities  as  in  Table  2. 


Pj-O 

.349(.oa) 

.13O(.00S) 

.0001. 004) 

.ieo(.oo4) 

.102  (.000) 

.240(.010) 

Pj-.S 

.36S(.017) 

.142 (.007) 

.070(.00l) 

.000(.004) 

.104 (.007) 

.205  (.017) 

Pj"** 

S 

1 

s 

.149(.009) 

.000(.004) 

,0001.004) 

.IIO(.OOO) 

.i00(.021) 

niKLI  3.  tstiSMitss  of  tbs  limitiiif  distributloii  for  ths  asmbor 
of  customers  in  qusos  tor  INS  euotomsrs  sa4 
flw  two  ssqusncos  of  sorwios  timM  sro  iaispsiitout 
prooMsos. 


Note  that  the  positive  correlation  of  both  sequences  of 
service  tiroes  tends  to  roake  the  probability  of  server  1 being 
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idle  and  the  probability  of  all  customers  being  at  server  1 
larger  than  ^en  only  one  sequence  of  service  tines  is 
positively  correlated.  The  change  in  the  limiting  distribution 
v(i)  is  again  somewhat  symnetric  in  i.  The  long  run  average 
number  of  custcmers  waiting  or  being  served  at  server  1 is 
2.5  if  @2  * ^2  * other  values  of 

( &2 ' ^2  ^ * 

It  is  clear  from  both  simulations  that  the  limiting 
distribution  for  number  of  customers  at  server  1 is  quite 
sensitive  to  serial  correlation  in  the  service  times.  The 
simulations  Indicate  that  perhaps  v(0)  « v(5)  in  the  case 
X ■ u « 1.  It  can  in  fact  be  shown  that  in  general 
Xtl  - v(0)l  » ytl  - v(N)]. 

5.  CONCLUSIONS 

In  this  paper  we  considered  one  scheme  for  using  EARMh 
processes  to  model  dependence  in  queues.  We  find  that  the 
introduction  of  dependence  does  affect  the  limiting  behavior 
of  the  queue.  There  are,  of  course,  many  other  schews  and 
Sana  of  these  will  be  oMsidesed  elseid»ere*  Two  advantages 
of  eeiag  the  IhMh  proeessss  ie  qeeuss  are  the  ease  of  intro- 
daelng  dspenfescs  in  the  qnees  and  the  ease  of  siamlating 
the  proossses.  The  maior  drawback  to  using  Bhmk  processes 

:i:  . 
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in  queues  is  the  difficulty  of  obtaining  the  exact  analytic 
results.  It  is  expected,  therefore,  that  approximation 
techniques  and  simulation  will  be  of  major  importance  in 
analysing  these  queues. 
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